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ABSTRACT 

Stars in galaxies form in giant molecular clouds that coalesce when the atomic hydrogen is converted 
into molecules. There are currently two dominant models for what property of the galactic disk 
determines its molecular fraction: either hydrostatic pressure driven by the gravity of gas and stars, 
or a combination of gas column density and metallicity. To assess the validity of these models, we 
compare theoretical predictions to the observed atomic gas content of low-metallicity dwarf galaxies 
with high stellar densities. The extreme conditions found in these systems are optimal to distinguish 
the two models, otherwise degenerate in nearby spirals. Locally, on scales < 100 pc, we find that the 
state of the interstellar medium is mostly sensitive to the gas column density and metallicity rather 
than hydrostatic pressure. On larger scales where the average stellar density is considerably lower, 
both pressure and shielding models reproduce the observations, even at low metallicity. We conclude 
that models based on gas and dust shielding more closely describe the process of molecular formation, 
especially at the high resolution that can be achieved in modern galaxy simulations or with future 
radio/millimeter arrays. 

Subject headings: Galaxies: dwarf - galaxies: ISM - ISM: molecules - stars: formation 



1. INTRODUCTION 

Theoretical arguments based on gravitational instabil- 
ity as well as observations of molecular gas reveal that 
low temperature (T ~ 10 K) and high density (n ~ 40 
cm~ 3 ) giant molecular clouds (GMCs) are the natu- 
ral sites where stars form. Although individual GMCs 
can be resolved only in the Milky Way or in a hand- 
ful of local galaxies (e.g. iBolatto et al.ll2008l and refer- 
ences therein), CO observations of several nearby spirals 
show that st ar formation mostly occurs in molecu l ar re- 
gionfl (e.g. IWong fc Blitz! [20021: iKennicutt etaLl [20071: 
Bigi el et al.ll2008D . At the same time, neutral atomic hy- 
drogen (H I) remains the primordial constituent of the 
molecular phase (H2), playing an essential role in the for- 
mation of new s tars, as shown by the low star-formation 
rate (SFR) (e. g. iBoselli fc G avazzi 2006) and low molec- 
ular content (|Fumagalli et al.l [200 9) found in H I-poor 
galaxies. Therefore, the transition from H I to H2 is a 
key process that drives and regulates star formation in 
galaxies. 

The problem of molecule formation has been stud- 
ied extensively in the literature mainly through two 
different approaches. The first is by modelling 
the formation of molecular gas empirically, start- 
ing from CO and H I maps in nearby galaxies. 
Followi ng this path, IWong fc Blitzl (|2002l hereafter 
WB02L Blitz fc Rosolowskvl (I2004L hereafter BR04), and 
Blitz fc Rosolowskvl ([2006L hereafter BR06) have inferred 
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that the molecular (£h 2 ) to atomic (Shi) surface density 
ratio 

Rh 2 = Sh 2 /Shi (1) 

in disks is a function solely of the hydrostatic midplane 
pressure P m , which is driven both by the stellar and 
gas density: i?n 2 ~ -Pm' 92 (hereafter BR model). The 
second approach models the microphysics that regulates 
the formation of H2 and its photodissociation. A de- 
tailed description should take into account the balance 
of H2 formation onto dust grains and its dissociation 
by Lyman- Werner (LW) photons and cosmic rays, to- 
gether with a complete network of chemical reactions 
that involves several molecules generally found in the in- 
terstellar medium (ISM). Due to this complexity, many 
studies address mainly the detailed physics of H2 in in- 
dividual clouds, without considering molecular forma- 
tion on galactic scales (e.g . the pioneering work by 
Ivan Dishoeck fc Bla,ckJll986l) . 

lElmegreenl (|1993l ) made an early attempt to produce 
a physically-motivated prescription for molecule forma- 
tion in galaxies by studying the H I to H2 transition 
in both self-gravitating and diffuse clouds as a function 
of the external ISM pressure P c and radiation field in- 
tensity j. This numerical calculation shows that the 



molecular fraction /h 2 = Sn 2 /5] ga 



P. 



with 



Sgas = Sh 2 +Shi- More recently, properties of the molec- 
ular ISM have been investigated with hydro dynamical 
simulations by iRobertson fc Kravtsovl (|2008f ) who have 
concluded that the H2 destruction by the interstellar ra- 
diation field drives the abundance of molecular hydro- 
gen and empirical relations such as the Rn 2 /P m corre- 
lation. Using numerical simul ations which i nclude self- 
consistent metal enrichme nt, IGnedin et al.l (|2009D and 
iGnedin fc Kravtsovl (|2010f ) have stressed also the impor- 
tance of metallicity in regulating the molecular fraction 
and therefore the SFR. Similarly. iPelupessv et al.l (2006) 
developed a subgrid model to track in hydrodynamical 
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simulations the formation of H2 on dust grains and its 
destruction by UV irradiation in the cold gas phase and 
collisions in the warm gas phase. 

A different approach based entirely on first 
principles has been pro p osed in a series of pa- 
pers bv IKrumholz et all ([IM hereafter KMT08), 
Krumholz et al.l (120091 hereafter KMT09), and 
McKee fc Krumhold (|2010l . hereafter MKfO). Their 



model (hereafter KMT model) describes the atomic- 
to-molecular transition in galaxies using a physically 
motivated prescription for the dust shielding and 
self-shielding of molecular hydrogen. This work differs 
from previous analyses mainly because it provides an 
analytic expression for /h 2 as a function of the total 
gas column density and metallicity (Z). Therefore, the 
KMT model can be used to approximate the molecular 
gas on galactic scales without a full radiative transfer 
calculation. 

In this paper we shall consider the BR and the KMT 
models as examples of the two different approaches used 
to describe the H I to H2 transition in galaxies. Re- 
markably, both formalisms predict values for /h 2 which 
are roughly consistent with atomi c and molecular obser- 
vations in local disk galaxies (see IKrumholz et al.l 120091 . 
sect. 4.1.3). The reason is that the BR model becomes 
dependent on the gas column density alone if the stellar 
density is fixed to typical values found in nearby galax- 
ies (see the discussion in Sect. [SJ. Despite the observed 
agreement, there are significant conceptual differences: 
the BR model is empirical and does not address the de- 
tails of the ISM physics, while the KMT model approxi- 
mates physically motivated prescriptions for the H2 for- 
mation as functions of observables. Hence, although in 
agreement for solar metallicity at resolutions above a few 
hundred parsecs, the two prescriptions may not be fully 
equivalent in different regimes. It is still an open question 
whether molecule formation is mainly driven by hydro- 
static pressure or UV radiation shielding over different 
spatial scales and over a large range of metallicities. 

A solution to this problem has important implications 
in several contexts. From a theoretical point of view, 
cosmological simulations of galaxy formation that span 
a large dynamic range will benefit from a simple prescrip- 
tion for molecular gas formation in order to avoid com- 
putationally intense radiative transfer calculations. Sim- 
ilarly, semi-analytic models or post-processing of dark- 
matter-only simulations will greatly benefit from a simple 
formalism that describes the molecular content in galax- 
ies. Observationally, the problem of understanding the 
gas molecular ratio has several connections with future 
radio or millimeter facilities (e.g. ALMA, the Atacama 
Large Millimeter Array or SKA, the Square Kilometer 
Array). In fact, these interferometers will allow high- 
resolution mapping of atomic and molecular gas across 
a large interval of redshift and galactic locations over 
which metallicity and intensity of the local photodissoci- 
ating UV radiation vary significantly. 

In this work, we explore the validity of the KMT 
and BR models in nearby dwarf starbursts both locally 
(< 100 pc) and on larger scales (~ I kpc). Their low 
metallicity (down to a few hundredths of solar values) 
combined with the relatively high stellar densities found 
in these systems offers an extreme environment in which 
the similarity between the two models breaks down. In 



fact, for a fixed gas column density, high stellar density 
corresponds to high pressure and therefore high molec- 
ular fraction in the BR model. Conversely, for a fixed 
gas column density, low metallicity in the KMT model 
results in a low molecular fraction (see Figured]). 

We emphasize that the BR model was not designed to 
describe the molecular fraction on scales smaller than 
severa l hundred parsecs; indeed, iBlitz fc Rosolowskvl 
(2006) explicitly warn against applying their model on 
scales smaller than about twice the pressure scale height 
of a galaxy. However, a number of theoretical models 
have extrapolated the BR model into this regime (e.g. 
iNaravanan et al.ll2010t: iMurante et al.l 12010ft . The anal- 
ysis we present at small scales (< 100 pc, comparable to 
the resolutions of simulations in which the BR model has 
been used) is aimed at highlighting the issues that arise 
from such an extrapolation. Furthermore, a compari- 
son of the pressure and shielding (KMT) models across 
a large range of physical scales offers additional insight 
into the physical processes responsible for the atomic-to- 
molecular transition. 

The paper is organized as follows: after a brief review 
of the two models in Section [3J we will present two data- 
sets collected from the literature in Section [31 The com- 
parison between models and observations is presented 
in Section [H while discussion and conclusions follow in 
Sections [5] and [6l Throughout this paper we assume a 
solar photosph e ric ab undance 12 + log(0/H) = 8.69 from 
lAsplund et all (|2009l ). Also, we make use of dimension- 



S gas /(M oP c- 2 ), stel- 



less gas surface densities E gas = 
lar densities /3 star = /9 s tar/(M0pc -3 ), gas velocity disper- 
sions i> gas = w gas /(km s _1 ), and metallicity Z' = Z/Zq. 

2. MODELS 

Here we summarize the basic concepts of the BR and 
KMT models which are relevant for our discussion. The 
reader should refer to the original works (WB02, BR04, 
BR06, KMT08, KMT09, and MK10) for a complete de- 
scription of the two formalisms. 

2.1. The BR model 

The ansatz at the basis of the BR model is that the 
molecular ratio i?H 2 is entirely determined by the mid- 
plane hydrostatic pressure P m according to the power- 
law relation: 



(2) 



The pressure can be evaluated with an approximate solu- 
tion for the hydrosta tic equilibrium fo r a two-component 
(gas and stars) disk (|Elmegreenlll989D : 



u gas 
^star 



(3) 



where £ s tar is the stellar density, and u s t ar and v gas are 
the stellar and gas velocity dispersion, respectively. 
For a virialized disk, equation (J3J reduces to 



P„ 



^frV 7, /y 0.5 i. -0.5 , y 0.5 fc-0.5\ . 
V (JrZjgagUgas I Lu Kas "- Kas \ '-'star star / ' 



(4) 



by replacing the surface densities with volume densities 
using the disk scale-heights h, the previous equation can 
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Fig. 1. — Left panel. Molecular fractions computed for the KMT (solid lines) and BR (dashed lines) models. Different lines represent 
three metallicities for the KMT model (from right to left, Z 1 = 0.1 blue, Z' = 1 black, and Z' = 10 red) or three stellar densities for the 



BR model (from right to left, p B 



- 0.001 blue, p s . 



: 0.1 black, and p 



10 red). For the KMT model, we assume a clumping factor 



c = 1. Blue compact dwarfs (BCDs) at low metallicities and high stellar densities are the optimal systems to disentangle between the two 
models which are degenerate in massive spiral galaxies with solar metallicity (compare the two black lines). Right panel. Models for the 
H I surface density as a function of the total gas column density, for the same parameters adopted in the left panel. While the KMT model 
exhibits a well defined saturation in the atomic hydrogen, in the BR model Shi increases asymptotically with E gas . (See the electronic 
edition of the Journal for a color version of this figure). 



be rewritten as 

Pm r 



„0.5 



gasi- 



Pstar; 



(5) 



Under the assumption that p s tar > Pgas, the contribution 
of the gas self-gravity can be neglected and equation ([5]) 
reduces to 



= 272 £„ 



(2p star )'°- 5 v 



K cm 3 



(6) 



where £ gas = S gas /(M Q pc 2 ), p : 
and 



/Ostar/(M pC 



-3\ 



gas 



;as /(km s _1 ). In deriving equation (|6]), con- 
stants in cgs units have been added to match equation 
(5) in BR06. Combining equation @ with equation (O, 
the molecular ratio in the BR model becomes 



i?H 2 



1.95 x 10" 



;as Pstar u gas 



'0.5 



V.,- 



0.92 



(7) 



where the best fit values P = (4.3 ± 0.6) x 10 4 K cm" 3 
and a = 0.92 have been derived from CO and H I obser- 
vations of local spiral galaxies (BR06). 

2.2. The KMT model 

The core of the KMT model is a set of two coupled 
integro-differential equations for the radiative transfer of 
LW radiation and the balance between the H2 phtodisso- 
ciation and its formation onto dust grains. Neglecting the 
H2 destruction by cosmic rays, the combined transfer- 
dissociation equation is 



V • F* = -na d cE* 



fmn 2 K 

/diss 



(8) 



On the left-hand side, F* is the photon number flux in- 
tegrated over the LW band. The first term on the right- 
hand side accounts for dust absorption, where n is the 
number density of hydrogen atoms, a d the dust cross 
section per hydrogen nucleus to LW-band photons, and 
E* the photon number density integrated over the LW 



band. The second term accounts for absorption due to 
photodissociation, expressed in term of the H2 forma- 
tion rate in steady state conditions. Here, /hi is the 
hydrogen atomic fraction while /di SS is the fraction of ab- 
sorbed radiation that produces dissociation rather than 
a de-excitation into a newly bound state. Finally, 1Z ex- 
presses the formation rate of molecular hydrogen on dust 
grains. 

Equation ([8|) can be integrated for a layer of dust and 
a core of molecular gas mixed with dust. This solution 
specifies the transition between a fully atomic layer and 
a fully molecular core and hence describes the molecular 
fraction in the system as a function of the optical depth at 
which F* = 0. Solutions to equation ([8]) can be rewritten 
as a function of two dimcnsionlcss constants 



and 



tr = na d R 



fdissO~dcE* 

nil 



(9) 



(10) 



Here, tr is the dust optical depth for a cloud of size R, 
hence equation Q specifies the dimensions of the system. 
Conversely, \ is the ratio of the rate at which the LW 
radiation is absorbed by dust grains to the rate at which 
it is absorbed by molecular hydrogen. 

To introduce these equations which govern the micro- 
physics of the H 2 formation into a formalis m that is appli- 
cable on galactic scales, one has to assume ()Wolfire et al.l 
2003) that a cold-neutral medium (CNM) is m pressure 
equilibrium with a warm- neutral medium (WNM). As- 
suming further that dust and metals in the gas com- 
ponent are proportional to the total metallicity (Z' = 
Z/Zq, in solar units), equations © and (fTU)) can be 
rewritten as a function of the observed metallicity and 
gas surface density. Using the improved formalism de- 
scribed in MK10, the analytic approximation for the 
molecular fraction as specified by the solutions of equa- 
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TABLE 1 

Data set for the high-resolution sample 81 



Name 


Distance 
(Mpc) 


12+log(0/H) 


— 

Mstar 

(10 5 M Q ) 


Radius 

(PC) 


- 

Pstar 

(M Q pc- 3 ) 




^star 

(M pc- 2 ) 


Shi 
(10 21 cm" 


S s f r 

~ 2 ) (MQyr-^pc" 2 ) 


IZwlS 


13 


7.19 


4.7 


56 


0.64 


47.9 


3.5 


0.134 








2.3 


56 


0.31 


23.3 






SBb0335-052 


53.7 


7.23 


10 


18.2 


39.5 


959 


7 


0.154 








11 




43.4 


1055 












A 
1 




15.8 


384 












11 




A Q A 

4o.4 


1055 












18 




71.0 


1726 












2 




7.9 


192 






Mrk71 


3.44 


7.90 


0.12 


7.2 


7.7 


74 


6 


0.214 








0.054 




3.5 


33 






UM462 


15.3 


7.98 


3.5-7.2 


21 


9-19 


255-520 


5.7 


0.060 








1.2-2.2 


27.5 


1.4-2.5 


51-92 












2.1 


43 


0.6 


36 












1.8 


27.5 


2 


74 












1.6-2.9 


33 


1-2 


48-85 












1.6 


50 


0.3 


21 






IIZw40 


10.3 


8.13 


6.4-12 


10.1 


149-280 


2005-3759 


7.9 


4.610 








1.3-15 


5.2 


222-2550 


1536-17730 






NGC5253 


3.5 


8.19 


0.13 


3.0 


115 


460 


6.1 


0.181 








0.7-4 


1.6-2.9 


3915-4080 


8704-15140 












10-13 


3.5 


5570-7240 


25980-33780 






NGC1140 


18.2 


8.20 


9.1 


7.3 


558 


5436 


2.5 


0.024 








59 


6.6 


4899 


43114 







a Individual references are provided in Appendix [A] together with a detailed description of how quantities are measured. 
Metallicities and SFRs are the same as those given in Table [2] 

b For the most uncertain values, we report the upper and lower limits. See Appendix [X] for further details. 



tion 



can be written as: 



/h 2 



1 



1 + 0.25s 



(11) 



for s < 2 and /h 2 = for s > 2. 
0.6 X + 0.01 X 2 )/(0.6r c ), X = 0.76(1 
Finally, E' 



= 0.066E comp Z'. 



comp 



Here s — ln(l + 
+ 3.1Z' 0365 ), and 
= cE gas where the 



clumping factor c > 1 is introduced to compensate for 
averaging observed gas surface densities over scales larger 
than the typical scale of the clumpy ISM. Primed surface 
densities are in units of M Q pc~ 2 . 

2.3. Differences between the two models 

In Figure Q] we compare the BR and KMT models 
to highlight some behaviours that are relevant to our 
analysis. In the left panel, we present molecular frac- 
tions computed using the KMT (solid lines) and BR 
(dashed lines) formalisms. Different lines reflect three 
choices of metallicity for the KMT model (from right 
to left, Z' = 0.1 blue, Z' = 1 black, and Z' = 10 
red) and three stellar densities for the BR model (from 
right to left, p stai = 0.001 blue, p stai = 0.1 black, and 
p star = 10 red). For a typical spiral disk with stellar 
mass M st ar = 10 10 M Q , size R — 10 kpc, and stellar 
height h = 300 pc, the stellar density is of the order 
of p stm . ~ 0.1. Figure [T] shows that, at solar metallicity 
and for a typical gas surface density E gas ~ 10—100, the 
two models predict similar molecular fractions. 

To break the degeneracy, we apply model predictions 
to observations of blue compact dwarf galaxies (BCDs) or 
low- metallicity dwarf irregulars (dlrrs) , characterized by 



high stellar density (p s 



1 — 100) and low metallicity 



(Z' = 0.3 — 0.03). In these environments, for a fixed gas 
surface density (excluding the limit /h 2 — > 1), the two 



models predict very different molecular to atomic ratios. 

In the right panel of Figure [TJ we show the predicted 
atomic gas surface density as a function of the total gas 
surface density, for the same parameters selected in the 
left panel. Besides the dependence on the metallicity and 
stellar density, this plot reveals a peculiar difference be- 
tween the two models. The KMT formalism exhibits a 
well defined saturation threshold in Ehi for a fixed value 
of metallicity. This corresponds to the maximum H I 
column density that is required to shield the molecular 
complex from the LW-band photons. All the atomic hy- 
drogen that exceeds this saturation level is converted into 
molecular gas. Conversely, the BR model has no satu- 
ration in the atomic gas surface density, but it increases 
slowly as the total gas surface density increases. 

3. THE DWARF GALAXY SAMPLES 

We study the behaviour of the KMT and the BR mod- 
els using two data sets compiled from the literature. 
Specifically, we have selected low-metallicity compact 
dwarf galaxies with sufficient observations to constrain 
gas densities, stellar masses, and metal abundances for 
a comparison with models. The first sample comprises 
16 BCDs and dlrrs, for which quantities integrated over 
the entire galaxy are available. These objects constitute 
a low-resolution sample with which we study the two 
models on galactic scales (> 1 kpc). For seven of these 
galaxies, we also have high-resolution H I maps and Hub- 
ble Space Telescope (HST) optical images. With these 
objects, we construct a high-resolution sample, useful to 
study the two formalisms at the scale of individual star 
cluster complexes (< 100 pc). 

Both models depend on the total gas surface den- 
sity. In principle, we could use CO emission, available 
in the form of integrated fluxes from the literature, to 
quantify E gas and the molecular content of individual 
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TABLE 2 

Data set for the low-resolution sample 



Name 


Distance 12+ 


log(0/H) a 


— 

Diameter 


M s t a r (min-avg) 




CO flux d 




SFR tracer 




(Mpc) 




(kpc) 


logM 


(M 0P c- 2 ) 


(K km/s) (M yr- 1 kpc- 2 ) 




Haro3 


16.8 


8.30 


5.23 - 4.79 


9.09 - 9.35 


10.2 


1.98 


0.106 


60fim 


IIZw40 


10.3 


8.12 


n on 1 n A 

2.80 - 1.04 


8.09 - 8.34 


91.4 


0.50 


4.610 


ft 


IZwl8 


13.0 


7.19 


1 A A J i ( i ■ ) 

1.44 - 0.93 


6.51 - 6.88 


39.3 


<1.00 


0.134 


ft 


Mrk209 


5.4 


7.81 


1.16 - 1.09 


7.29 - 7.36 


14.2 


0.45 


0.051 


Ha 


Mrk33 


24.9 


8.45 


7.61 - 6.86 


9.59 - 9.79 


4.3 


6.21 


0.103 


60/xm 


Mrk71 


3.4 


7.90 


0.81 - 0.81 


6.54 - 6.89 


12.4 


<0.34 


0.214 


Ha 


NGC1140 


18.2 


8.20 


6.56 - 6.55 


9.30 - 9.58 


15.7 


0.97 


0.024 


Ha 


NGC1156 


7.1 


8.23 


5.53 - 5.93 


8.62 - 9.10 


5.0 


0.76 


0.007 


Ha 


NGC1741 


55.1 


8.05 


17.82 -22.44 


9.29 - 9.81 


9.7 


1.53 


0.040 


60/xm 


NGC2537 


8.0 


8.19 


4.48 - 3.72 


9.14 - 9.26 


5.4 


0.56 


0.013 


Ha 


NGC4214 


3.3 


8.20 


2.91 - 7.19 


8.66 - 8.67 


2.3 


0.90 


0.003 


Ha 


NGC5253 


3.5 


8.19 


3.29 - 3.14 


8.63 - 9.05 


2.7 


0.73 


0.181 


ff 


NGC7077 


13.3 


8.04 


2.72 - 2.90 


8.54 - 8.60 


3.6 


0.68 


0.014 


Ho 


SBS0335-052 


53.7 


7.23 


4.24 - 3.35 


7.79 - 8.72 


24.5 


<5.43 


0.154 


Ha 


UM448 


81.2 


8.00 


15.67 - 8.18 


10.30 - 10.69 


34.1 


0.82 


0.669 


60/xm 


UM462 


15.3 


7.97 


2.48 - 2.44 


7.95 - 8.20 


16.1 


0.55 


0.060 


60/xm 


a References 


for metallicity: 


Davidgc 


(1989); Campos- Aeuilar et al. 


(1993); llzotov et al.l (1997); 


Gil dc Paz et al] (120031: 


Izotov & Thuan (2004a): Thuan 


& Izotov (2005): van Zee & Havncs (2006). 










b Diameters as 


computed from stellar profil 


es and derived from NED. 











c References for H I: HyperLeda except for IZwl8 and Mrk71, as described in the text. 
d Refe r ences for CO:ISage et alj Jl99^:ILeon et alj tl99ll:ITavlqr et"al] (fl99l ) ; iBarone et all (200Cft : lDale et al.l (pOOll ') : IGil de Paz et al.l 
pOu^l t lATbrecht et al. l 120041): ILerov et al.) 1200511; ILerov et al.l 120071). 

e References for SFR: lYmmg et alj (119891): IDrissen et al.l 1200 0): Hop kins et al.l 12002T ) ; IGil de Paz et all (200% ; IHunt et all j2004); 
llzotov fc Thuanl (2004a); H unt et al.l (2005a); H unt et ali i 2005b): Schmit t*etaIT $006): van Zee & Havncs (2006). 



galaxies. However, recent studies of mole cular hydrogen 
trace d through a gas-to -dust ratio (e.g. Ilmara fc Blitd 
I2007t ILerov et al.l I2009T ) support the idea that CO is 
a poor tracer of molecular hydrogen in low metallicity 
environments, most ly due to its inability to self-shield 
(jWolfire et al.ll201Ch . Therefore, CO seems an unreliable 
H2 tracer for these metal-poor galaxies. For this unfortu- 
nate reason, we avoid any attempt to precisely quantify 
Sh 2 j but rather use the observed H I column density as 
a lower limit on the total gas column densities. 

As discussed in Section l2~5I the KMT model has a well- 
defined saturation threshold for Ehi, and this threshold 
constitutes an observationally-testable prediction. The 
BR model does not have such a threshold and at a given 
Pstar is in principle capable of producing arbitrarily high 
values of Shi provided that the total gas density E gas 
is sufficiently high. However, the extremely weak varia- 
tion of Sri with E gas at large total gas column density 
(Shi <x ^gas 8 ) means that the amount of total gas re- 
quired to produce a given Ehi may be implausibly large. 
This effect allows us to check the BR model as well us- 
ing only H I (Section fJJ, albeit not as rigorously as we 
can test the KMT model. We also check the robustness 
of our results in Appendix [C] where we impose an up- 
per limit on Eh 2 eit her from SFRs, as suming a depletion 
time £dc P i ~ 2 Gyr (Bi giel et al.ll2008t ) for molecular gas, 
or from CO fluxes, using a conservative CO-to-H2 con- 
version. 

In the next sections, we discuss in detail the proce- 
dures adopted to derive gas surface densities and stellar 
densities for the two samples. The reader not interested 
in these rather technical aspects can find the analysis, 
discussion and conclusions starting from Section 2J 

3.1. High-resolution sample 

Seven BCDs are found in the literature with high- 
resolution H I maps and with sufficient ancillary HST 



data to infer stellar masses on scales < 100 pc, typical of 
individual GMCs. A detailed description of how we com- 
pute Ehi and p s tar in individual galaxies is provided in 
Appendix [A] together with a list of relevant references. 
Here, we only summarize the general procedures we use. 

Stellar masses of individual clusters are in a few cases 
directly taken from the literature. Otherwise, we infer 
stellar masses from integrated light by comparing two 
methods. The first is based on age estimates, whenever 
those are available in the literature. In this case, we in- 
fer stellar masses from observed absolute magnitudes by 
comparing the K or V band luminosity with predictions 
at th e given age by Starburst99 (SB99; iLeitherer et al.l 
1999). This is done assuming an instantaneous burst, 
similar metallicity, and a Salpeter initial mass function 
(IMF) with lower and upper mass limits at 1 and 100 M Q , 
respectively. The second method is based on optical and 
near-infrared colors (e.g., B — V , V — I, V — J, V — H, 
V — K). In this case, we use mass-to- light ( M/L) ra- 
tios inferred from colors (|Bell &; de JonelfeOOlD . and the 
stellar masses are derived directly from observed lumi- 
nosities. Usually, the two methods give similar results to 
within a factor of ~ 2. 

Once the masses are known, we obtain stellar densities 
with sizes taken from the literature. If not available, 
we measure them by fitting a two-dimensional elliptical 
Gaussian to the clusters in HST images. For the closest 
objects, in order to avoid resolving individual stars, we fit 
binned surface brightness profiles with a one-dimensional 
Gaussian. 

Stellar masses are probably the most uncertain quan- 
tities in our study. In fact, our first method suffers from 
the rapid changes in the broadband output of a starburst 
at young ages (4 — 10 Myr), due to the onset of red su- 
pergiants whose amplitude and time of onset depends on 
metallicity Moreover, the initial mass function of the 
SB99 models and the lower-mass cutoff may introduce 
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TABLE 3 

Photometric quantities for the low resolution sample 3. 



Name B K K ([3.6]) b K ([4.5]) c B — K (min.) B — K (mean) B — K (std) M star , K (min.) M starjK (avg) 

(mag) (mag) (mag) (mag) (mag) (mag) (mag) (M ) (M©) 



Haro3 


13. 


22 


10.61 




10.05 


2.61 


2 


.89 





.40 


9.09 


9.35 


IIZw40 


11, 


91 


10.98 


10.89 


10.22 


0.96 


1 


.24 


0. 


. 12 


8.09 


8.34 


IZwl8 


15 


.86 


15.24 


14.59 


14.45 


0.62 


1 


.10 





.42 


6.51 


6.88 


Mrk209 


15. 


.09 


12.50 




12.68 


2.40 


2 


.19 





.12 


7.29 


7.36 


Mrk33 


13. 


.39 


10.42 


10.09 


10.13 


2.96 


3 


.17 





.18 


9.59 


9.79 


Mrk71 


11. 


.44 




11.74 


10.90 


-0.30 





.12 





.60 


6.54 


6.89 


NGC1140 


13. 


46 


10.48 


9.92 


10.07 


2.98 


3 


.30 





.29 


9.30 


9.58 


NGC1156 


11. 


.35 


9.45 


8.49 


8.51 


1.90 


2 


.53 





.55 


8.62 


9.10 


NGC1741 


13 


08 


11.82 


10.76 


10.78 


1.26 


1 


.96 





.61 


9.29 


9.81 


NGC2537 


12 


.09 


9.11 


8.91 


9.02 


2.98 


3 


.07 





.10 


9.14 


9.26 


NGC4214 


10. 


.15 


7.90 


7.92 




2.22 


2 


.23 





.02 


8.66 


8.67 


NGC5253 


10. 


.63 


8.21 


7.60 


7.23 


2.41 


2 


.95 





.50 


8.63 


9.05 


NGC7077 


13. 


.90 


11.36 


11.39 


11.47 


2.42 


2 


.19 





.06 


8.54 


8.60 


SBS0335-052 


16. 


.45 


15.40 


14.02 


12.79 


1.05 


2 


.38 


1 


.31 


7.79 


8.72 


UM448 


14. 


,11 


11.33 


10.65 


10.49 


3.11 


3 


.62 





,15 


10.30 


10.69 


UM462 


14. 


.42 


12.70 


12.33 


12.20 


1.72 


2 


.01 





.26 


7.95 


8.20 



density limit found by [Hopkins et al.l (|2010j ) , and most 
are 5 — 10 times below this limit. Interestingly, the stellar 
densities here are uncorrelated with metallicity, implying 
that some other parameter must play the main role in 
defining the properties of massive star clusters. 

3.2. Low-resolution sample 

We have collected a second sample from the literature 
by requiring only that quantities integrated over the en- 
tire galaxy be available. Due to the lower spatial resolu- 
tion, this data set is suitable to study the KMT and BR 
models on larger scales (> 1 kpc). Our search yielded a 
total of 16 low-metallicity star-forming galaxies; among 
these are the 7 objects in the high-resolution sample. 
We have compiled gas and stellar densities, distances, 
and metallicity for these 16 objects, most of which are 
classified as BCDs, but some are dlrrs (Sm, Im), since 
they are more diffuse, larger in size, and more luminous 
(massive) than typical BCDs. 

Stellar masses are compu ted from Svitzer /IRAC fluxes 
following the formulation of lLee et al.l (|2006fl , as we sum- 
marize in Appendix IBl together with a comment on the 
dominant sources of uncertainty. Stellar densities are 
then derived assuming spherical symmetry and the sizes 
inferred from the stellar component, as measured from 
IRAC images. The resolution of these images (~1"2) 
is a factor of 10 lower than the worst HST resolution, 
so that the compact regions are unresolved. This implies 
that the stellar densities derived for this sample are much 
lower than the values quoted for individual star clus- 
ter complexes. Moreover, for non-spherical (spheroidal) 
BCDs these densities correspond formally to lower lim- 
its; the volume of a prolate spheroid is smaller than the 
volume of a sphere by a factor (6/a) ' 5 , with a and b 
the semi-major and semi-minor axis, respectively. In our 
sample, the mean axis-ratio is a/b ~ 1.5 with a 0.5 stan- 
dard deviation. This discrepancy is small enough to jus- 
tify our assumption of spherical symmetry. In any case, 
the possible volume overestimate could partially compen- 
sate the potential overestimate of stellar density because 
of nebular emission contamination or free-free emission 
(see the discussion in Appendix [Bj . 



a See Apeendix [Bj for a detailed description of the listed quantities. 
b Assumed K - [3.6] = 0.03 
c Assumed K - [4.5] = 0.20 

additional uncertainty, up to a factor o f 3, considering a 
full range of systematic uncertainties (jBell et al.l f2003). 
Instead, sources of error in the second method are the 
strong contribution of nebular continuum and line emis- 
sion to the broadband colors of young starbursts. This 
can be a particularly severe problem in the K band be- 
cause of recombination lines and free-free emission which 
in some cases constitutes as much as 50% of the broad- 
band K m agnitude (see iVanzi etUhl 120001 : iHunt et all 
I2001L [2003h . Despite this rather large uncertainty on the 
stellar densities, the results presented in the next sections 
can be considered rather robust. In fact, a variation in 
the density larger than the uncertainty would be required 
to significantly alter our conclusions. A more extensive 
discussion on this issue is presented in Section 14.1.21 

To complete the data set, we add to the gas and stellar 
densities values for the metallicity, distances, and SFR 
indicators as collected from the literature. Individual 
references are provided in Appendix [X] and Tables [1] and 
[21 We derive integrated SFRs using Ha, 60/im and radio 
free-free fluxes as different tracers. The fina l rates are 
given assuming the empirical calibrations by Kennicutt 
(|1998ft for the Ha, bylHopkins et al.l (|2002D for the 60Atm, 
and by IHunt et ah! (|2005al) for the radio free-free emis- 
sion. We note that this last tracer is optimum in the ab- 
sence of non-thermal emission, as typical in young star- 
bursts. Since SFRs are used to set an upper limit on 
the total gas density assuming a given depletion time 
(Appendix [Cj , we choose the maximum value whenever 
more than one indicator is found for a single galaxy. To- 
tal SFR surface densities are then calculated adopting 
the galaxy sizes from NEE0- A summary of the collected 
and derived data is presented in Table [TJ 

The stellar densities in the high-resolution sample 
are generally quite high, and associated with massive 
compact star clusters, some of which are in the Su- 
per Star Cluster (SSC) category (e.g.. lO'Connell et al.l 
fl99l : lMeurer et alJll995HWhitmore e t al. 2005]). Despite 
their extreme properties, none of the BCDs in the high- 
resolution sample exceeds the maximum stellar surface 

4 NASA/IPAC Extragalactic Database. 
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Fig. 2. — Comparison between observations of the H I surface 
density (crosses) with prediction from the KMT model (lines) as 
a function of the metallicity. Different curves are computed for 
different total gas column densities (£ gas = 50 solid line, S gas = 
10 2 dotted line, E gas = 5 X 10 2 dashed line, and £ gaa = 10 3 dash- 
dotted line). The KMT model is computed for c = 1 and H I 
data are shown as lower limits since they are averaged over regions 
which are more extended than 100 pc (see the text for further 
details). Observations do not rule out the KMT model since there 
is a significant overlap between the parameter space allowed by 
both the theoretical curves and the data. 

For most object s, integrated H I fl uxes are retrieved 
from HyperLedcf] (jPaturel et al. 2003). We then convert 
integrated fluxes into mean column densities using opti- 
cal radii from NED, and assuming that the gas extends 
twice as far as the stellar component (se e ILee et aLll2002t 
Ivan Zee et al.lll998atlThuan et al.|[200l . For I Zw 18, the 
integrated H I flux is not availab le in HyperLeda and 
we co nsider the flux published in Ide Vaucouleurs et al.l 
(1991;). Similarly, for Mrk71 we estimate the total 
atomic gas from available inter ferometric observa tions 
averaged over the entire galaxy (Thuan "eTaII[200l . 

12 CO(l — 0) fluxes are available for most of the galax- 
ies here considered (see Table [2]). For three galaxies, the 
most metal-poor objects in our sample (SBS 0335— 052 E, 
IZwl8, and Mrk71), we find only CO upper limits 
in the literature. Because we use CO fluxes only to 
set upper limits on Eh 2 (Appendix [Cj) . we choose one 
of the largest CO -to-H2 conversion measured to date 
(jLeroy et a l. 2009) . It is worth noting that for extremely 
metal poor galaxies (e.g. 1 Zw 18) the adopted conversion 
factor may still underestimate the H2 content. To make 
our limits even more conservative, we compare these val- 
ues with Eh 2 inferred from SFRs and we choose for each 
galaxy the maximum of the two. 

As with the high-resolution sample, we derive SFRs 
from Ha, 60/im, or free-free emission (see Table[2]for ref- 
erences). Again, SFR densities are computed assuming 
the optical size as given by NED. Finally, we collect infor- 
mation on the metallicity and distances for each object. 
A summary of the data derived for the low-resolution 
sample is given in Table [2] 

4. ANALYSIS 
5 http://leda.univ-lyonl.fr 



4.1. Testing models on small scales (< 100 pc) 

With the aim of testing how the BR and KMT models 
perform at high stellar density and low metallicity, we 
first compare both formalisms with the observed H I sur- 
face densities and stellar densities in the high-resolution 
data set. Since the average resolution of this sample is 
below 100 pc, this part of the analysis focuses mainly on 
the molecular fraction in individual GMCs and associa- 
tions rather than on larger ISM spatial scales. As pre- 
viously mentioned, the BR formalism was not developed 
to describe the molecular fraction in such small regions 
(BR06). Hence the results presented in this section are 
intended to assess possible pitfalls of extrapolating the 
pressure model to small scales. Furthermore, a compar- 
ison of the performances of the BR and KMT models 
below 100 pc provides insight into what quantities are 
relevant to the production of molecules on different size 
scales. 

As summarized in Section [2~2l the KMT formalism de- 
scribes the molecular fraction as a function of the total 
gas column density and metallicity. A free parameter is 
the clumping factor c that maps the observed column 
density £ gas onto the relevant quantity in the model, 
i.e. the column density of the cold phase in individual 
clouds. With resolutions coarser than ^100 pc, beam 
smearing dilutes the density peaks and one must adopt 
c > 1 in order to recover the intrinsic gas surface density 
Scomp = cS gas . However, given the high resolution of 
the HST images, we set c = 1 so that the KMT model 
has no free parameters. 

Conversely, as reviewed in Section 12.11 the BR model 
describes i?n 2 as function of the total column density and 
stellar volume density. An additional parameter in this 
case is the gas velocity dispersion, set to v gas = 8. Apart 
from a similar dependence on £ gas , a direct comparison 
between models and observations is not straightforward. 
We start our analysis by confronting each model with 
observations, and then attempt a comparison of both 
models and data. 

4.1.1. The KMT model predictions below 100 pc 

In Figure [2] we present the observed H I surface den- 
sity (crosses) together with predictions from the KMT 
model (lines) as a function of the metallicity. Differ- 
ent curves are computed for different total gas column 
densities (£ gas = 50 solid line, £ gas = 10 2 dotted line, 
^gas — 5 x 10 2 dashed line, and E gas = 10 3 dash-dotted 
line). Here, and for the rest of this analysis, we correct 
the gas column density for helium with a standard coef- 
ficient 1.36. We do not include corrections for projection 
effects because in dwarf galaxies a unique inclination an- 
gle is not well defined for a warped (non-planar) H I 
distribution or in a triaxial system. 

As discussed in Appendix \K\ interferometric H I obser- 
vations do not achieve the resolution required to match 
the HST observations. A possible solution is to down- 
grade HST images to match the atomic hydrogen maps. 
However, since the exact value for c would be unknown at 
the resultant resolution (> 100 pc), we perform our anal- 
ysis on scales < 100 pc, compatible with HST images and 
where c — > 1. For this reason, we express the observed 
H I column density as lower limits on the local Shi- This 
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is because coarser spatial resolutions most likely average 
fluxes on larger areas, thus lowering the inferred peak col- 
umn density. Indeed, whenever H I observations at differ- 
ent resolutions are compared, better resolution is associ- 
ated with higher column densities. SBS 0335—052 E is an 
example: this BCD has Eht = 7.4 x 10 20 cm" 2 in a beam 
of 2075 X 157 (IPustilnik et aljUOOll) and 2x 10 21 cur 2 in 
a 3"4 beam (|Ekta et al.l 120091 ). The inferred H I column 
density is even higher, 7 x 10 21 cm -2 , with the smaller 
2" beam in HST/GHR S observations of Lya absorption 
(|Thuan fc IzotovlH997l ). 

In any case, the lower limits illustrated in Figure [2] 
prevent us from concluding that model and observations 
are in complete agreement, although this is strongly sug- 
gested. Adopting a conservative approach, this compar- 
ison shows that observations do not immediately rule 
out the KMT model on scales of < 100 pc; 5/7 of the 
galaxies here considered are consistent with predicted 
curves. Although not crucial for the current and re- 
maining analysis, the quoted metallicity may in some 
cases overestimate the dust and metal content which 
contributes to the H2 formation. In fact, in the KMT 
model, it is the CNM that plays a relevant role in reg- 
ulating /h 2 an d the assumption that the nebular metal- 
licity reflects the metal abundances in the cold ISM 
may not hold in all cases. Specifically, the optically- 
inferred metallicity used here is dominated by the ion- 
ized phase. Studies of the metal enrichment of the neu- 
tral gas in metal poor dwarfs show that the neutral 
phase can be some times less metal-enriched than the ion- 
ized medium (e.g. iThuan et al.ll2005t iLebouteiller et al.1 
120041: iLecavelier des Etangs et al.l 120041 ) . " Furthermore, 
galaxies with the lowest nebular metallicities have sim- 
ilar neutral gas abundances, while dwarfs with higher 
ionized nebular metallic ities can have up to ~7 times 
(|Lebouteiller et aT1l2009D lower neutral ISM abundances 
(see however iBowen et al.ll2005t ) . A recent interpretation 
for this effect is that although mixing is effective in diffus- 
ing new metals from ionized reg ions, due to the larger vo l- 
ume, the enrichment is modest (|Lebouteiller et aTj 2009). 
This justifies the use of a single metallicity for multiple 
GMCs since, in the worst case, we would overestimate 
by some factor the local metal content. Data points, es- 
pecially the ones at higher metallicity, would be offset to 
lower values and the parameter space common to data 
and model would increase, mitigating the discrepancy 
found at Z' 0.5. 

4.1.2. The BR model predictions below 100 pc 

Turning our attention to the BR model, we test predic- 
tions (lines) against the observed Shi (crosses) in Figure 
[H Stellar densities are available for individual GMCs 
measured from high-resolution HST images. Therefore, 
while multiple observations overlap in metallicity in Fig- 
ure [2] here we show distinct data-points for the same 
galaxy. Once again, different curves are for a selection of 
total gas density (E gas = 50 solid line, E gas = 10 2 dotted 
line, E' . = 5x 10 2 dashed line, = 10 3 dash-dotted 

line, and E gas = 10 4 long-dashed line); observed H I 
column densities are represented as lower limits on the 
local atomic gas column density. When stellar densities 
are particularly uncertain (see discussion in Appendix 
[A"| . we plot both lower and upper limits connected with 
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p„„(M s pc- 3 ) 

Fig. 3. — Comparison between the predicted Ejji from the BR 
model (lines) and observations (crosses) as a function of the stellar 
density in individual GMCs and associations. Different curves are 
for choices of total gas density (S gaa = 50 solid line, S gas = 10 2 
dotted line, E_.„ = 5 X 10 2 dashed line, £„„ = 10 3 dash-dotted 

line, and S gas = 10 4 long-dashed line). Observed H I column 
densities are represented as lower limits on the local atomic gas 
column density (see the text for further details). Observations rule 
out the BR model once extrapolated below 100 pc. 

a dotted line. 

According to the BR model, despite the low metallicity, 
a high fraction of hydrogen is expected to be molecular 
because of the enhanced stellar density. However, Figure 
[3] illustrates that observations discourage the use of pres- 
sure models on scales < 100 pc. Even under the very 
conservative hypothesis that stellar densities are over- 
estimated by a factor of 2 — 3 and that the total col- 
umn densities can reach very high values (e.g. the long- 
dashed curve at E gas = 10 4 ), observations mostly lie in 
the region not allowed by the extrapolation of the BR 
modejj. Compatibility between the extrapolation of the 
BR model and a good fraction of the data would require 
Egas 10 10 M pc -2 . A value that large would corre- 
spond to Ay > 10 5 with a dust-to-gas ratio that is 1% 
of the Milky Way value, and is thus ruled out by the fact 
that the star clusters are observable. Moreover, such a 
large E gas would make the gas mass in the observed re- 
gion larger than either the baryonic or the dark matter 
mass of the entire dwarf galaxy. Clearly, even though wc 
cannot directly detect the molecular component of the 
gas, we can rule out the presence of such a large amount 
of gas on other grounds, and we can therefore conclude 
that the extrapolated BR model is incompatible with the 
observations. We give a more rigorous estimate of the 
maximum plausible value of E gas in Appendix C. 

An additional tunable parameter in the BR model is 
the gas velocity dispersion and a substantial change in 
v gas can affect its predictions. In this paper, following 

6 One might attempt to improve the agreement data-model by 
artificially smearing the stellar density down to the same resolution 
as the H I observations (while ignoring the presence of additional 
stars outside the HST PSF). This brings the total gas surface den- 
sity required to match the majority of the observations down to 
S ~ 10 3 — 10 4 , still a very large value. 
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BR04 and BR06, we adopt v gas = 8. Since pressure varies 
linearly with the velocity dispersion, we can solve for the 
value of v gas required for the BR model to match the ob- 
servations. We find that < 2 in order to have one half 

gas ~ 

of the data points consistent with the mo del; this is in 
contrast with recent H I observations (e.g. iChung et al.l 
120091 iWalter et al.ll2008l ) that show typical dispersion ve- 
locities v gas > 5 (and in many cases v gas > 10) in all the 
surveyed galaxies. It is worth mentioning that the use 
of the observed v gas is not always appropriate, although 
unavoidable; for example, w gas depends on the thermal 
velocity and the gas in a cold medium has a lower ve- 
locity dispersion than what inferred from a multiphase 
ISM. We conclude that the disagreement found in Figure 
[3] cannot be explained with uncertainties on the velocity 
dispersion. 

Finally, we should assess if the high pressure predicted 
by the model can be attributed to the use of hydrostatic 
equilibrium in a disk rather than in a sphere, which would 
be more appropriate for our systems. Intuitively, this is 
not the case since the central pressure in a sphere of gas 
and stars cannot be lower than the midplane pressure 
of the disk. In fact, the central point in a sphere has 
to support the weight of the entire system, while each 
point in the midplane of a disk has to support only the 
pressure from the components along the vertical direc- 
tion. This argument is substantiated by a quantitative 
analysis. A solution of the hydrostatic equilibrium equa- 
tion for a sphere of gas and stars shows that the cen- 
tral gas pressure is enhanced by a quantity that depends 
on (i>star/^gas) 2 - Therefore, to minimize an increase in 
the pressure due to the stellar component, the condition 
(^star/vgas) ~ 1 nas to be satisfied. However, an increase 
in v gas is reflected by an increase in the gas pressure itself. 
In other words, the central pressure in a gas sphere with 
cold kinematics (low v gas ) receives a significant contribu- 
tion from the stellar potential {v s t ar /v gas > 1), while a 
gas sphere with hot kinematics (high v gas ) has an intrin- 
sically higher gas pressure (vstar/vgas < 1). 

4.1.3. A direct comparison between the two models 

A different way to visualize both the models and the 
observations for the high-resolution sample is shown in 
Figure 21 where we plot predictions for Shi of the KMT 
model (solid lines) and the BR model (dashed lines) as a 
function of the total gas column density. Different curves 
in the KMT model are for the maximum, central, and 
minimum metallicity observed in the sample (from the 
bottom to the top, Z' = 0.32 red, Z' = 0.12 black, and 
Z' = 0.03 blue). The curves for the BR model corre- 
spond to the maximum, central, and minimum stellar 
density (from the bottom to the top, \ogp stm . = 3.86 

red, log^tar = 1-38 black, and log/o star = —0.52 blue). 
As in the previous figures, observed lower limits on Shi 
are superimposed. The yellow-shaded region in Figure |4] 
indicates the maximum parameter space allowed by the 
observations. This plot summarizes the two main results 
presented in the previous paragraphs. Observations of 
Shi reveal that an extrapolation of the BR model below 
scales of 100 pc results in a significant overestimation of 
the molecular fraction. In fact, for the observed Shi, ex- 
ceedingly high total gas surface densities (S > 10 4 ) 
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Fig. 4. — KMT model (solid lines) and BR model (dashed lines) 
as a function of the total gas density. Different curves are for the 
maximum, central, and minimum observed values of metallicity (in 
the KMT model: Z' = 0.32 red, Z' = 0.12 black, and Z' = 0.03 
blue) and stellar density (in the BR model: logp star = 3.86 red, 

logp star = 1.38 black, and logp star = —0.52 blue). Observed lower 
limits on Sjji are superimposed. The yellow shaded region encloses 
the largest parameter-space allowed by observations. Data rule out 
any extrapolation of the BR model below 100 pc, where it largely 
overestimates the molecular fraction. Conversely, observations do 
not rule out immediately the KMT model. (See the electronic 
edition of the Journal for a color version of this figure). 
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Fig. 5. — Comparison of the models with data in the low- 
resolution sample. Solid lines are for the KMT formalism, for the 
maximum (Z 1 = 0.58 red) central (Z' = 0.20 black) and mini- 
mum [Z' = 0.03 blue) observed metallicity. Dashed lines are for 

the BR model, for the maximum (logp star = —1.45 red), central 

(logp star = —2.00 black), and minimum (logpstar = —3.18 blue) 
stellar density in the sample. Lower limits on the total column den- 
sity are for £ ga s = Shi- At low resolution, stellar densities drop 
by orders of magnitude and both models appear to be consistent 
with the data. (See the electronic edition of the Journal for a color 
version of this figure). 

are required by the BR model to reproduce observations. 
As shown quantitatively in AppendixO such high values 
appear to be unrealistic. Conversely, observations seem 
to suggest a good agreement between the KMT model 
and data. Also, comparing Figure [4] with Figure [IJ it 
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appears that the different behaviour of the two models 
for a similar gas column density is related to which quan- 
tity regulates the molecular fraction at the second order, 
subordinately to the gas column density. In fact, the 
discrepancy with observations and the BR model is as- 
sociated with the high values of stellar densities, while 
the consistency between the observed H I column den- 
sities and the KMT model is fostered by the low value 
of metallicity that raises the atomic hydrogen saturation 
limit. 

4.2. Testing models on galactic scales (> 1 kpc) 

In the second part of this analysis, we compare predic- 
tions from models and observations on larger scales (> 1 
kpc), by considering spatially integrated quantities for a 
larger sample of BCDs. Before we start, it is worth men- 

> 20 (see BR04) which 
3 holds also for the 



tioning that the condition S s 
ensures the validity of equation 
low- resolution data set. 



4.2.1. A comparison between models and global data 

In Figure [5] we present a comparison between observed 
H I surface densities and models, as previously done in 
Figure [4] for the high-resolution sample. Solid lines rep- 
resent the KMT model, for the maximum (Z 1 = 0.58 
red), central (Z 1 = 0.20 black), and minimum (Z 1 = 0.03 
blue) observed metallicity. Dashed lines are for the BR 
model for the maximum (log/9 star = —1.45 red), cen- 
tral (log/9 star = —2.00 black), and minimum (logp star = 
—3.18 blue) stellar density. Lower limits on the total gas 
column density are computed for S gas = Shi (see Ap- 
pendix [C] for a version of this figure that includes upper 
limits). 

A comparison of Figures [5] and [4] reveals that the BR 
model predicts for the galaxy as a whole a much higher 
H I surface density, compared with the predictions for 
regions smaller than 100 pc. By going from the high- 
resolution sample to the low-resolution one, we lose the 
ability to analyse the local structure of the ISM, and 
are limited to average quantities which dilute the den- 
sity contrasts in both gas and stars over many GMC 
complexes and aggregations. As a result, Shi and p s t ar 
are lowered by one and two orders of magnitude respec- 
tively. This behaviour is reflected in the BR model as a 
decrease in the pressure by an order of magnitude, which 
now guarantees an overall agreement between data and 
the model. Conversely, the beam-smearing is accounted 
for in the KMT model by the clumping factor, here as- 
sumed to be c ~ 5. Hence, despite the different spatial 
scales, the KMT formalism is able to account for the 
mean observed Shi- 

4.2.2. A test with individual galaxies 

Although indicative of a general trend, Figure [5] does 
not allow a comparison of individual galaxies with mod- 
els. This is particularly relevant for the KMT formal- 
ism that predicts a saturation in the H I column density 
as a function of the metallicity. To gain additional in- 
sight, we present in Figure |6] the predicted H I surface 
density (Shi, mod) against the observed one (Shi,ods) for 
individual galaxies. For the KMT model, we compute 
Shi, mod using observed metallicities (red crosses). Since 
the total gas column density is unknown, we compute 




10' 

W (M pc 2 ) 

Fig. 6. — Predicted H I surface densities (Shi, mod) computed ei- 
ther with observed metallicity using the KMT model (red crosses) 
or with stellar density using the BR model (blue squares) are shown 
as a function of the observed values (Em,obs) in individual galax- 
ies. For the KMT model, lines connect lower and upper limits on 
S H j mo< j, derived assuming proper limits on S gas from the satura- 
tion threshold. The BR model is consistent with observations in 
individual galaxies, even at low metallicity. Similarly, the KMT 
formalism can account for most of the galaxies, but fail to repro- 
duce observations in some cases (See the electronic edition of the 
Journal for a color version of this figure). 




Fig. 7. — Ratio of the observed and predicted H I surface den- 
sities in the KMT model (red crosses) as a function of metallicity 
(upper panel) and stellar density (lower panel). The two dashed 
horizontal lines indicate the perfect agreement between data and 
model. There are no evident systematic trends in the discrepancies 
between data and the KMT formalism (See the electronic edition 
of the Journal for a color version of this figure). 



for each galaxy a range of Shi, mod (shown with solid red 
lines) using upper and lower limits on the observed total 
gas surface density. Lower limits are derived assuming 
Sgas = Shi,odsj while upper limits arise self-consistently 
from the H I saturation column density, naturally pro- 
vided by the model. For the BR model, we instead com- 
pute only lower limits on Shi, mod with stellar densities 



(blue squares), assuming S g 



Sm,obs (upper limits 



derived for Shi, mod are presented in Appendix |C|) . 
From Figure [6l the asymptotic behaviour of Shi in 
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Fig. 8. — Comparison between the observed H I surface density 
and the value predicted by the KMT model, as in Figure [6] but 
for an arbitrarily lower metallicity (Z' = Z' /3). A possible cause 
for the observed discrepancy data-model is a lower metallicity in 
the cold phase than the one observed in the ionized gas. (See the 
electronic edition of the Journal for a color version of this figure). 

the BR model prevents a tighter constraint on £ gas . In 
general, there is good agreement between observations 
and model predictions for all the galaxies, despite the low 
mean metallicity of this sample. Conversely, the KMT 
model allows a narrower interval of H I surface density 
because of the well defined atomic hydrogen saturation. 
Therefore Figure [6] provides a more severe test of the 
KMT formalism that nevertheless reproduces correctly 
most of the observations. Although for 4/16 galaxies the 
KMT predictions are inconsistent with the data, there 
seems to be no peculiar reason for the failure of the model 
for these objects (see Section |4.2.3|) . 

For many objects the agreement between models and 
observations occurs close to the lower limits on Eni.mod) 
i.e. when E gas = Ehi- We stress that this is not an 
obvious outcome of the assumption made on the total 
gas surface density, since high stellar densities in the BR 
model or high metallicity in the KMT model would imply 
a high molecular fraction irrespectively of Em.obs- If we 
change perspective for a moment and assume that mod- 
els are a reliable description of the molecular hydrogen 
content, the observed trend suggests that low-metallicity 
galaxies are indeed H I rich. Perhaps, this is not surpris- 
ing due to the reduced dust content at low metallicity 
and the consequent reduction of the shielding of molec- 
ular gas from the LW-band photons. 

4.2.3. Systematic effects for the KMT model 

Finally, in Figure [7] we explore whether the KMT 
model exhibits systematic effects within the range of val- 
ues allowed by the low-resolution sample. For this pur- 
pose, we present the ratio EHi,obs/Effl,mod f° r * ne KMT 
model (red crosses) as a function of the observed metal- 
licity (top panel) and stellar density (bottom panel). As 
in Figure El we display with solid lines the full interval of 
Eui.mod- The lack of any evident trend either with metal- 
licity or stellar density suggests that the KMT model is 
free from systematic effects. In particular, because the 
four deviant galaxies are not found at systematically high 
or low metallicity, the difficulties in assessing metallicity 



for the cold-phase gas could be the cause of such de- 
viations. This hypothesis was touched upon in Section 
14.1.11 where we comment on the possibility that metals 
in the cold gas can in some cases be several factors lower 
than that observed in the ionized gas. An example of 
the importance of a correct metallicity determination is 
provided in Figure [SJ where we repeat the comparison 
between the KMT model and observations in individual 
galaxies after we arbitrarily redefine Z' = Z' /3. As ex- 
pected, for lower metallicity, the H I saturation moves 
to higher atomic surface densities and the KMT model 
better reproduces the observations. 

5. DISCUSSION 

The analysis presented in Sect. 0] indicates that ex- 
trapolations of the BR model based on pressure overpre- 
dict the molecular fraction on small spatial scales (< 100 
pc), while this formalism recovers the observed values on 
larger ones. Such a failure of the BR model on small 
scales was predicted by BR06, who point out that the 
model does not properly account for the effects of gas self- 
gravity or local variations in the UV radiation field. Con- 
versely, the KMT model based on gas and dust shielding 
is consistent with most of the observations both locally 
and on galactic scales, although it is more prone to ob- 
servational uncertainties in the ISM structure (through 
the clumping factor) and the cold-phase metallicity. In 
this section, after a few comments on these points, we 
will focus on an additional result which emerges from 
our comparison: the molecular fraction in galaxies de- 
pends on the gas column density and metallicity, while 
it does not respond to local variations in pressure from 
enhancements in the stellar density. 

5.1. The effect of self- gravity at small scales 

The problem of gas-self gravity is that it introduces an 
additional contribution to the force balance on scales typ- 
ical for GMCs. In this case, pressure equilibrium with the 
external ISM is no longer a requirement for local stability 
and the empirical power law in equation ([7]) may break 
down. However, since self-gravity enhances the internal 
pressure compared to the ambient pressure in equation 
(|5|). one would expect locally even higher molecular frac- 
tions than those predicted by an extrapolation of the BR 
model. This goes in the opposite direction of our results, 
since the observed molecular fraction is already overesti- 
mated. Hence, a different explanation must be invoked 
for the data-model discrepancy in Figure 

5.2. The effect of the radiation field 

A reason for the high Ehi observed on small spatial 
scales is related to the intensity of the UV radiation field. 
In fact, regions which actively form stars probably have 
an enhanced UV radiation field compared to the mean 
galactic value. Since the BR formalism does not explic- 
itly contain a dependence on the UV radiation field in- 
tensity, it is reasonable to expect discrepancies with ob- 
servations. In contrast, the KMT model attempts to 
explicitly account for local variation in j, by considering 
how such variations affect conditions in the atomic ISM. 
This makes it more flexible than the BR model in scal- 
ing to environments where conditions vary greatly from 
those averaged over the entire galaxy. 
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We stress here that the BR model is not completely 
independent of the radiation field, but simply does not 
account for a variati on in j. Assumin g the scaling rela- 
tion /h 2 ~ p 2 2 j~ 1 (jElmegreenl 119931 ). the BR formal- 
ism is commonly considered valid only when variations 
in pressure are much greater than those in j, allowing 
us to neglect the latter. Being empirically based, the 
BR model contains information on a mean j, common 
for nearby spirals. Therefore this model describes the 
molecular content as if it were only regulated by pres- 
sur^H- Hence, it can be applied to describe the molecular 
fraction only on scales large enough such that variations 
in the local UV intensity are averaged over many com- 
plexes and, eventually, j approaches a mean macroscopic 
value similar to that found in the galaxies used to fit the 
BR model. 

This idea is quantitatively supported by recent numer- 
ical simulati ons. The molecular fraction o f the galaxies 
simulated by iRobertson fc Kravtsovl (|2008f) is consistent 
with the observed Rn 2 ~ P 92 only when the effects 
of the UV radiation field are taken into account. In 
fact, when neglecting the radiation field, a much shal- 
lower dependence Rn 2 ~ P 0A is found. In their dis- 
cussion, the observed power-law index a ~ 0.9 results 
from the combined effects of the hydrostatic pressure 
and the radiatio n field. Starting from /h 2 oc p 2 2 j~ 1 
(|Elmegreenlll993D . assuming a Kennicutt-Schmidt law in 
which j ~ Sgf,. ~ £ gas , and under the hypothesis that the 
stellar surface density is related to the gas surface density 
via a star formation efficiency £ s tar ~j^L S ; -^h 2 = P a 
requires that ( Roberts on fc Kravtsov|[2008l ) 



a = 2.2 



(1 + /?/2) - 
1 + 0/2 



(12) 



In the simulated galaxies, different star formation laws 
and efficiencies (mostly dependent on the galaxy mass) 
conspire to reproduce indexes close to the observed a ~ 
0.92, in support of the idea that a mean value for j is 
implicitly included in the empirical fit at the basis of the 
BR model. 

5.3. The effect of stellar density 

So far, we have discussed why the BR model extrapo- 
lated to small spatial scales is unable to predict the ob- 
served H I surface density due to its reliance on a fixed 
"typical" j. Galaxies in our sample are selected to be 
metal poor, but some of them have a higher SFR (median 
~ 0.6 M© yr _1 ) and higher specific star formation rate 
(SSFR; median ~ 10 -9 yr _1 ) than observed in nearby 
spiralfS (< 10- 10 yr" 1 ; e.g. IBothwelllTal1l200lfl . If the 
UV intensity were the only quantity responsible for the 
disagreement between observations and the BR model 
below 100 pc, we would also expect some discrepancies 
on larger scales for galaxies with enhanced star forma- 
tion. However, such discrepancies are not observed. As 

7 Incidentally, this exact statement can be found in BR06. 

8 We recall our choice of the maximum SFR among the values 
available in the literature to conservatively obtain an upper limit 
on the molecular emission. For this reason, our median SFR is 
biased towards high values. In any case, a SFR ~ 0.3 Mq yr — x , 
typical of BCDs (Hopkins ct al. 2002), would correspond to a SSFR 
5 x 10~ 10 yr" 1 at the median mass M sta r = 6 x 10 8 M Q . 




Fig. 9. — Atomic gas surface density predicted by the BR model 
(black lines) as a function of the stellar density. Different curves 

are for different values of the total gas surface density (£ gas = 10 
solid line, £ gas = 15 dotted line, S gaa = 50 dashed line, 2 gas = 10 2 
dash-dotted line, and S gas = 10 3 long-dashed line). The two hori- 
zontal lines are for the KMT model at the mean metallicity in the 
high-resolution sample (red-dotted line) and in the low-resolution 
sample (blue-dashed line). The open diamond and the open trian- 
gle represent the median values of stellar density and H I surface 
density in the high-resolution and low-resolution samples, respec- 
tively. The open square represents a typical stellar density and H I 
surface density for local spiral galaxies, from a median of the values 
in table 1 of BR04, assuming a stellar disk height h = 300 pc. Our 
analysis favours a model in which the molecular fraction in the ISM 
depends to second order on metallicity rather than stellar density. 
(See the electronic edition of the Journal for a color version of this 
figure) . 

discussed, on larger scales, both the KMT and the BR 
models are able to reproduce observations, despite the 
different assumptions behind their predictions. 

We argue that there is an additional reason for the ob- 
served discrepancy between data and the extrapolation 
of the BR model below 100 pc. While both models de- 
pend to first order on the gas column density, our analysis 
favours a model in which the local molecular fraction in 
the ISM depends to second order on metallicity rather 
than density as in the BR model. To illustrate the ar- 
guments in support of this hypothesis, in Figure [9] we 
show the atomic gas surface density predicted by the BR 
model (black lines) as a function of the stellar density. 
Different curves are for different values of the total gas 
surface density (S gas = 10 solid line, E gas = 15 dotted 
line, £ gas = 50 dashed line, E gas = 10 2 dash-dotted line, 
and S gas = 10 3 long-dashed line). Superimposed, there 
are three data points. The open diamond and the open 
triangle represent the median values of stellar density and 
H I surface density for the high-resolution (S HI = 65, 

Pltar = 16 ) and low-resolution (E HI = 15, p' stax = 0.01) 
samples, respectively. The open square represents a typ- 
ical stellar density and H I surface density (S HI = 7, 

Pstar = 0-2) for local spiral galaxies, from a median of 
the values in table 1 of BR04, assuming a stellar disk 
height h = 300 pc. Finally, the two horizontal lines are 
for the KMT model at the mean metallicity in the high- 
resolution sample (red-dotted line, Z — 0.12) and in the 
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low-resolution sample (blue-dashed line, Z = 0.20). 

Recalling that for the BR model P ~ £ ga s%isv / Pstan 
we see from Figure [9] that the stellar density does not 
provide a major contribution to the variation in pres- 
sure when p s tar ^ 0.5, typical for large galactic regions 
(open triangle and square). In this regime, the predicted 
Ehi from the BR model (black lines) is mostly depen- 
dent on varaition in the total gas column density alone. 
In fact, for a constant velocity dispersion, the pressure 
model becomes equivalent to models based on gas shield- 
ing. This reconciles on theoretical grounds the equiva- 
lence observed in local spirals between the BR and KMT 
(blue horizontal dashed line) models. At local scales, 
i.e. moving towards higher stellar density, p s t ar provides 
an important contribution to the total pressure and the 
extrapolation of the BR model predicts high molecular 
gas fractions, as seen in Figure HO The expected Ehi 
drops accordingly, following a trend that is in disagree- 
ment with observations (open diamond) . Conversely, the 
KMT model is insensitive to p s t a r and accounts only for 
a variation of the gas column density and metallicity. In 
this case, the predicted Ehi (red horizontal dotted line) 
moves towards higher values in agreement with observa- 
tions. 

From this behaviour, we conclude that stellar density 
is not a relevant quantity in determining the local molec- 
ular fraction. Furthermore, because the KMT formal- 
ism recovers the high observed H I column density for 
low metallicity, dust and metals have to be important in 
shaping the molecular content of the ISM . This is further 
supported bv lGnedin fe Kravtsovl (|2010[) who show with 
simulations how the observed star formation rate, which 
mostly reflects the molecular gas fraction, depends on 
the metallicitj0. It follows that the observed dependence 
on the midplane pressure is only an empirical manifesta- 
tion of the physics which actually regulates the molecular 
fraction, i.e. the effects of the UV radiation field and the 
gas and dust shielding. 

5.4. Fixed stellar density for molecular transitions 

A final consideration regards the observational evi- 
dence that the transition from atomic to molecular gas 
occurs at a fixed stellar density with a small variance 
among different galaxies (BR04). This result favors 
hydrostatic pressure models because if the atomic-to- 
molecular transition were independent of stellar surface 
density, there would be much more scatter in the stel- 
lar surface density than is observed. However, this 
empirical relation can also be qualitatively explained 
within a formalism based on UV radiation shielding. In 
fact, if the bulk of the star formation takes place in 
molecule dominated regions, the build-up of the stellar 
disks eventually will follow the molecular gas distribu- 
tion, either directly or via a star formation efficiency (se e 
iRobertson fc Kravtsovll2008t iGnedin fc Kravtsovll2010D . 
Since the transition from molecular to atomic hydro- 
gen occurs at a somew hat well-defined ga s column den- 
sity (KMT09; WB02; iBigiel et al.l 12003 ), it is plausi- 
ble to expect a constant surface stellar density at the 

9 However, a detai l ed ana lysis of star forming regions in BCDs 
by IHira shita & Hunt ( 2004) suggests that other parameters such 
as gas density, size, and geometry play a role in determining the 
local SFR. 



transition radius. While this picture would not apply 
to a scenario in which galaxies grow via subsequent 
(dry) mergers, rece nt hydrodynamical simulations (e.g. 
iBrooks et al.l [2009) support a model in which stars in 
disks form from in-situ star formation from smoothly ac- 
creted cold or shock-heated gas. 

6. SUMMARY AND CONCLUSION 

With the aim of understanding whether the princi- 
pal factor that regulates the formation of molecular 
gas in galaxies is the midplane hydrostatic pressure 
or shielding from UV radiation by gas and dust, we 
compared a pressure model (BR mo del; iWong fc Blitz! 
[200alBlitz fc Rosolowskvl[2^M[2006l ) a nd a model based 
on UV pho t odissociation (KMT mode l; iKrumholz et al.1 
120081 120091 : iMcKee fc Krumholzl |20ToT ) against observa- 
tions of atomic hydrogen and stellar density in nearby 
metal-poor dwarf galaxies. Due to their low metallic- 
ity and high stellar densities, these galaxies are suitable 
to disentangle the two models, otherwise degenerate in 
local spirals because of their proportionality on the gas 
column density. 

Our principal findings can be summarized as it follows. 

- On spatial scales below 100 pc, we find that an 
extrapolation of the BR model (formally applica- 
ble above ~ 400 pc) significantly underpredicts the 
observed atomic gas column densities. Conversely, 
observations do not disfavour predictions from the 
KMT model, which correctly reproduces the high 
H I gas surface densities commonly found at low 
metallicities. 

- Over larger spatial scales, with the observed and 
predicted H I surface density integrated over the 
entire galaxy, we find that both models are able to 
reproduce observations. 

- Combining our results with numerical simula- 
tions of the molecular formation in the galax- 
ies IS M (jElmegreenl 119891 : IRobertson fc Kravtsovl 
2008) which indicates how the UV radiation field 
(j) plays an essential role in shaping the molecu- 
lar fraction, we infer that the discrepancy between 
the BR model and observations on local scales is 
due partially to the model's implicit reliance on an 
average j, which breaks down at small scales. In 
contrast, the KMT model properly handles this ef- 
fect. 

- Since on scales ~ 1 kpc the BR model agrees with 
observations despite the low metallicity and high 
specific SFR in our sample, we infer that the dis- 
crepancy between pressure models and observa- 
tions below 100 pc also arises from their depen- 
dence on stellar density. An increase in stellar den- 
sity corresponds to an increase in the hydrostatic 
pressure which should, in the BR model, reduce the 
atomic gas fraction. No such trend is seen in the 
observations. 

- If we drop the dependence on the stellar density, 
the pressure model reduces to a function of the to- 
tal gas column density and becomes equivalent to 
the KMT model, for a fixed velocity dispersion and 
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metallicity. This provides a theoretical explanation 
for the observed agreement of the two models in lo- 
cal spirals. 



In conclusion, our analysis supports the idea that the 
local molecular fraction is determined by the amount of 
dust and gas which can shield H2 from the UV radia- 
tion in the Lyman- Werner band. Pressure models are 
only an empirical manifestation of the ISM properties, 
with the stellar density not directly related to the H2 
formation. Although they are useful tools to character- 
ize the molecular fraction on large scales, obviating the 
problem of determining the clumpy structure of the ISM 
or the metallicity in the cold gas as required by models 
based on shielding from UV radiation, pressure models 
should be applied carefully in environments that differ 
from the ones used in their derivation. These limitations 
become relevant in simulations and semi-analytic mod- 
els, especially to describe high-redshift galaxies. Fur- 
thermore, a correct understanding of the physical pro- 
cesses in the ISM is crucial for the interpretation of ob- 
servations, an aspect that will become particularly rele- 
vant once upcoming facilities such as ALMA will produce 
high-resolution maps of the ISM at high redshifts. 

Combining our analysis with both theoretical and ob- 
servational efforts aimed at the description of the ISM 
characteristics and the SFR in galaxies, what emerges 
is a picture in which macroscopic (hence on galactic 
scales) properties are regulated by microphysical pro- 
cesses. Specifically, the physics that controls the atomic 
to molecular transition regulates (and is regulated by) 



the SFR, which sets the UV radiation field intensity. The 
ongoing star formation is then responsible for increasing 
the ISM metallicity and building new stars, reducing and 
polluting at the same time the primordial gas content. 
Without considering violent processes more common in 
the early universe or in clusters, this chain of events can 
be responsible for a self-regulated gas consumption and 
the formation of stellar populations. Future and ongoing 
surveys of galaxies with low-metallicity, active star for- 
mation and high gas fraction (e.g., LITTLE THINGS; 
Hunter et al.) will soon provide multifrequency obser- 
vations suitable to test in more detail the progress that 
has been made on a theoretical basis to understand the 
process of star formation in galaxies. 
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APPENDIX 
NOTES ON INDIVIDUAL GALAXIES 
IZw 18 

The main body of I Zw 18 consists of two main clusters, the north-west (NW) and the south-east (SE) components 
with an angular separation of ~ 6". A third system, known as "Zwicky's fiare" or "C" component, lies about 22" to 
the northwest of the NW cluster. H I maps are available from Ivan Zee et aLl (|1998a|) . together with the HST/WFPC2 
F814W image (0'.'045 resolution). The H I peaks close to t he fainter SE cluste r, rather than to the NW where the stellar 
density is higher. Lya observations with HST/GHRS bv lKunth et al.l ([1994) (2"x 2" beam) are also available for the 
NW c loud and have a better resolution than the VLA map. At the assumed distance of 13 Mpc (|Izotov fc Thuanl 
2004b), 1" = 63 pc. Stellar masses of the two massive clusters in I Zw 18 are not pu blished. Hence, multiband 
integrated photometry of the two star clusters in IZw 18 is taken from lHunt et al.l ((2003), together with cluster ages 
as modeled by them. Sizes of the clusters are measured from fitting Gaussians to the surface brightness profiles (not 
previously published). 

The SE cluster has an age of 10 Myr and, near where the distribution peaks (see Ivan Zee et al.l[l998ah . M K = -12.4 
(in a 2" aperture). The lowest-metallicity SB99 models (Z' = 0.001) give Mk = —12.4, which implies a stellar mass 
of 2.3 x 10 5 M Q . With a K-band lumino sity of 1.91 x 10 6 L Q , this would give a (M/L) K = 0.12, as inferred from the 
SB99 models. The lBell fc de Jongl (pOOlT ) predictions give (M/L) K = 0.17, on the basis ofV-K and (M/L) K = 0.09, 
from V — J. Hence, the SB99 value of 0.12 is roughly consistent. We therefore adopt the value of 2.3 x 1O 5 M for 
the stellar mass of the SE cluster . We can check th e inferred stellar mass by inspecting the if-band surface brightness 
at the SE peak (see Figure 5 in iHunt et al.l l200"3l ) . /zk = 19.2 mag arcsec -2 . This gives £k = 185.3 L pc~ 2 and, 
assuming (M/L)k = 0.12 from SB99, we would have S star = 22.2 M pc" 2 . This is in good agreement with the value 
of S s tar = 23.3 M pc -2 , inferred from the absolute luminosity (see above) and the measured radius of 56 pc. 

The NW cluster has an age of 3 Myr and Mk = —13.25 (in a 2" aperture). The lowest-metallicity SB99 models give 
Mk = —16.2, which is rather uncertain because of the rapid increase in luminosity at about 3 Myr w hen the most 
mass ive stars start evolving off of the main sequence, a phase which is not correctly described in models ■: Oriulia et al.1 
1999). In fact, the observed V — H color of 0.29 is predicted by SB99 to occur at ~10 Myr, not at 3 Myr, which is the 
best-fit photometric age. In any case, the inferred mass from this model is 6.6 x 1O 4 M . With a iC-band luminosity 
of 1.91 x 10 6 L , this would give a (M/L)k = 0.016. The same exercise repeated for the V band, with My = —12.86 
and the SB99 prediction of M v = -15.7, give an inferred mass of 7.3 x 10 4 M o , and (M/Z) v = 0.006. These M/L 
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values are quite low, roughly 6 times smaller than those predicted bv lBell fc de J ong (2001) from the observed colors 
of I Zw 18. We therefore use the latter M/L ratio. With V - K = 0.38, V — H = 0.29, and V — I= -0.04, we estimate 
(M/L)k = 0.11, 0.10, and 0.09, respectively. Therefore, adopting 0.10, we derive a stellar mass of 4.2 x 1O 5 M . 
Repeating the calculations for V band, we find (M/L)y = 0.039, 0.034, and 0.033, respectively. Adopting 0.033, 
we would derive a similar stellar mass of 3.9 x 10 5 M Q . Again, the if-band surface brightness of the NW cluster 
(^ K = 18.3 mag arcsec -2 , 0'.'5 resolution) gives a similar result. We find Ek = 536 L Q pc -2 , and, with (M/L) K = 0.10, 
becomes 53.6 M pc~ 2 . With a cluster radius of 56 pc (0'.'89), this would correspond to a cluster mass of 5.3 x 10 5 M , 
about 1.3 times that inferred from the lower-resolution photometry. Hence, to obviate problems of resolution (1" radius 
aperture or 63 pc, vs. a 56 pc radius measured from the HST image), we adopt the mean of these two measurements 
for the stellar mass of the NW cluster, namely 4.7 x 10 5 M Q . 

SBS 0335-052 E 

SBS 0335— 052 E hosts six SSCs, with most of the star format ion activity centered on the two brightest ones to the 
southeast. The H I distribution is pub lished in lEkta et all (|2009D . and the HST/ ACS F555M image (0"050 resolution), 
was published by Rei nes et al.l (|2008| ). The H I map is of relatively low resolution (~ 3'.'4) and does not resolve the 
six SSCs individual l y sinc e they are distributed (end-to-end) over roughly 2'.'6. Lya observations with HST/GHRS 
bv IThuan fe Izotovl (Il997h (2"x2" beam) are also available. In our analysis, we use this column density, being at 
better resolution than the one derived from H I emission map. At the assume d distance of 53.7 Mpc, 1" = 260.3 pc. 
Stellar masses for individual clusters have been derived by Reines et al. (2008|) by fitting the optical and UV spectral 
energy distributions, and we adopt these masses here. Comparison with masses inferred from if-band is unfruitful 
since nebular and ionized gas contamination make this estimate highly uncertain. We measure the size of the clusters 
by fitting two-dimensional Gaussians. They are unresolved at the HST/ ACS resolution of 0'.'050, but since they have 
the same size to within 13%, we assume the average radius of 18.2 pc. Therefore, the inferred mass densities result in 
lower limits. 

MRK11 

Mrk71 (NGC 2363) is a c omplex of H II regio ns in a larger irregular galaxy, NGC 2366. There are two main knots of 
star- formation act i vity ( see IDrissen et al.ll2000D . denoted A and B. A low-resolution H I map (12'.'5x ll'/5) is available 
from lThuan et all (|2004j ). but at this resolution we are unable to distinguish the tw o main clusters which are 5" apart. 
At the assumed distance of 3.44 Mpc (derived from Cepheids. iTolstov et al.lll995l) . 1" = 16. 7 pc. Stellar masses of the 
two star burs t knots in Mrk 71 are no t published. Hence, I^-band photometry wa s taken from | Drissen "eTalllpOOOT ). and 
/-band from IThuan fc Izotovl (|2005l) , together with cluster ages as modeled bv IDrissen et al.l ( 20001 ) . As for IZwl8, 
sizes of the clusters were measured from fitting ID Gaussians to the surface brightness profiles. 

The knot A has V = 17.3 mag and, after correcting for Ay = 0.3 mag, we derive an absolute magnitude My — —10.4. 
At an age of 3 Myr, SB99 models (at Z=0.00 4) predict My = -15.2 We would thus infer a stellar mass of 1.2 x 10 4 M , 
and an implied (M/L)y ratio of 0.012. ThelBell fe de Jond (|2001ft predictions give (M/L)y = 0.010, on the basis of 
stella r V - K = -0.42 as modeled bv lNoeske et all (|2000f l and {M/L) v = 0.10, from stellar B — V = -0.19 (also 
as in INoeske et all [20001) . The latter value from B — V is a factor of 10 higher than the former from V — K , and 
highly inconsistent with the SB99 value of 0.01. The /-band photometry of knot A from lThuan fe Izotovl (|2005l ) gives 
a similar inconsistency. With / = 17.97 mag, and a corre sponding absolute magnitude of 9.71, we would infer a stellar 
luminosity of 3.43 x 10 5 L . With lBell fe de Jond pOOl (M/L)i values of 0.025 (from V-K) and 0.15 (from B- I/), 
we would derive stellar masses of 8.6 x 10 3 and 5.1 x 10 4 M Q , respectively. Since three values are roughly consistent 
(~ 10 4 M Q ), we adopt 1.2 x 10 4 M Q as the mass for knot A. 

Knot B is slightly older than knot A (4 Myr) consisten t with its Wolf-Rayet stars and strong stellar winds as inferred 
from P Cygni-like profiles in the UV (Dri ssen et al.ll2000l) . It is also slightly fainter with V = 18.05 (after correcting for 
Ay = 0.3 mag), corresponding to an absolute magnitude My = —9.63. SB99 models (at 4 Myr) predict My — —15.3, 
which would give a stellar mass of 5.4 x 1Q 3 M , and an implied {M/L)y ratio of 0.005. Again, we derive M/L 
ratios as a function of color from iBell fc de J ong (2001), and obtain {M/L)y — 0.06 from stellar V — K = 0.67 and 
(M/L)y = 0.15, from stellar B — V = -0.07 (|Noeske et alll2000D . With a stellar V-band luminosity of 6.1 x 10 5 L Q , 
we would infer a stellar mass of 3.6 x 10 4 M Q with (M/L)y = 0.06, and 9.2 x 10 4 M o w ith (M/L)y = 0T5. Both 
masses are larger than those inferred for knot A, inconsistently with the observation of IDrissen et al.l (|2000f) that 
knot B contains only ~6% of th e ionizing photons n ecessary to power the entire H II region. Nevertheless, knot A is 
suppo sedly enshrouded in d ust (|Drissen et alll2000l) . so the situation is unclear. The /-band photometry of knot B 
from lThuan fc Izotovl ([2005) is not edifying. With 7" = 18.92 mag, and a corresponding absolute magnitude of —8.76, 
we would infer a stellar luminosity of 1.43 x 10 5 L Q . With lBell fc de Jond pOOll ) (M/L)i values of 0.10 (from V-K) 
and 0.20 (from B — V), we would derive stellar masses of 1.4 x 1O 4 M and 2.9 x 1O 4 M , respectively. These values 
are all greater t han the mass inferr ed for knot A, even though knot B is reputed to be intrinsically ~16 times fainter 
(see above, and IDrissen et al1l2000l ). For this reason, we adopt the SB99 value of 5.4 x 10 3 M as the mass for knot B. 

UM462 

U M462 hosts six SSCs, with 9" separation from end to end. The H I distribution is available in the literature 
from Ivan Zee et al.l (|1998b[) . together with the ground-based ESO/SOFI Ks image (0'.'28 pixels, 0'.'8-l"seeing) from 
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iVanzil (|2003f) . The resolution of the H I map (6'/6x5'.'2) is just barely sufficient to distinguish the two clusters. At 
the assumed distance of 13.5 Mpc, 1" = 65.4 pc. Stellar masses of the six SSCs in UM 462 have been derived by 
IVanzil (|2003l ) by estimating the age from the Ha equivalent width, then comparing SB99 models at that age with the 
Ha luminosity of the cluster after correcting for extinction. These estimates differ from the other SB99 comparisons 
described for previous galaxies because they extend to a lower lower -mass l imit, 0.1 M Q rather than 1 M . We thus 
consider a range of possible masses given by the values published bv IVanzil ([20031 ) and what we infer from comparing 
the -Ks-band luminosities of the individual clusters with S B99 p redictions as above (at the published age). On the 
basis of the obs e rved V - Ks and V - J colors (jVanzi et alj|200a IVanzi 2003]), we d erive (M/L)k = 0.05 according to 
IBell fe de Jorigl (12001 . A comparison of these numbers with the values provided by IVanzil (120031) reveals t hat ty pical 
uncertainty on the mass calculations are roughly a factor of 2 or less. Sizes are measured bv IVanzi" etal] ([2002T I. but 
the clusters are unresolved at the ground-based resolution. Hence, the mass densities are formally lower limits. 

IIZw40 

IIZw40 is a cometary BCD with two tails. The main star-formation activity is occurring at the "head of the 
co met" , namely in tw o knots in a north-south orientation, separated vertically by 1/5. The upper knot, dubbed "A" 
bv IVanzi et al.l (|2008l ), is elongate d alon g roughly an east- west direction and contains the rising-spectrum thermal 
radio sources found by Beck et al. (2002). The lower B knot i s round, fainter than A, and apparently does not host 
any compact radio sources. The H I distribution is published in Ivan Zee et al.1 (|1998bD and HST/F814W images (0'.'025 
pixels) are also available. The resolution of the H I map, 5'.'7x4'.'8, does not distinguish the two knots seen at HST 
resolution. At the assumed distance of 10.3 Mpc, 1" = 49.9 pc. IIZw40 is located sufficiently near the plane of the 
Milky Way that the foreground extinct ion is quite high, Ay = 2.7 mag. Stellar masses of the two clusters in IIZw40 
have been derived bv IVanzi et al.l (|2008l ). but with a Kroupa IMF and using Br7 luminosity (see also UM462) in a 15" 
aperture. Moreover, a subtraction procedure was applied that made assumptions about the Br7 flux of region A in the 
vicinity of knot B. Hence, to test these values for stellar masses, we recompute them from continuum measurements, 
using similar procedures to those us ed for the preyious galaxies. 

Ages for the clusters are taken from lVanzi et al.l ()2008l ) , who compared observed Br7 equivalent widths to SB99 model 
predictions. Photometry in the HST F555W, F814W, and F160W passbands is performed by fitting a 2D Gaussian 
to each knot. This determination is roughly consis tent with aperture photometry in a 1" aperture. Correcting the 
observed magnitudes for the high extinction as in IVanzi et al.l ()2008f ) from the NIR hydrogen recombination lines 
(Ay = 4.0 mag for knot A, Ay = 4.9 mag for knot B), we have 14.62, 15.80, 13.34 for F555W, F814W, and F160W, 
respectively (knot A), and 14.13, 14.58, and 13.31 for F555W, F814W, and F160W (knot B). Under the approximation 
F555W~V , we derive an absolute magnitude of —15.4 for knot A, to be compared with the SB99 prediction of —15.2 
(at 3 Myr), and —15.9 for knot B to be compared with —15.5 (at 7 Myr). The comparisons with SB99 give stellar 
masses of 1.2 x 10 6 M Q and 1.5 x 1O 6 M , for knots A and B, respectively (the inferred (M/L)y are 0.01 and 0.007 
for knots A and B). With the approximation F555W~V , F8 14W~I, and F160W~ H, V -I = -1.18 and -0.45, and 
V- H = 1.31 and 0.82 for knots A and B, respectively. Using IBell fe de Jongl (|2001[ L we derive (M/L)y = 0.0005 and 
0.004 from V — I, and (M/L) v = 0.21 and 0.09 from V — H. 

These values differ substantially from those inferred from SB99. In particular, the value of 0.0005 is unrealistic, and 
difficult to reconcile with other galaxies and other M/L inferred for IIZw40. Hence, we consider a range of stellar 
masses for knot A with 6.4 x 10 5 M Q (for (M / L)y = 0.005) and 1.2 x 1O 6 M (SB99) as lower and upper limits. For 
knot B instead, we assume a lower limit at 8.1 x 10 5 M (for (M/L)y = 0.004) and an upper limit at 1.5 x 10 6 M 
(SB99). These are roughly consistent with the masses given bv IVanzi etaLl (|2008ft of 1.7 x 1O 6 M for knot A and 
1.3 x 10 5 M for knot B, as inferred from comparing SB99 predictions of Br7 emission over a 15" aperture (knot A) 
and 0"75 aperture (knot B). For the sizes of the two clusters, we fit a ID Gaussians to the surface-brightness profiles in 
the F814W band and obtain 10.1 pc for knot A and 5.2 pc for knot B. The dimensions derived by fi t ting 2 D Gaussians 
are smaller, namely 4.4 pc (A) and 3.7 pc (B), similar to the dimensions obtained bv iVanzfe t al. (200j|) for the star 
cluster itself (rather than the more extended H II region emission) . 

NGC5253 

NGC 5253 is a nearby dwarf galaxy in the Centaurus group at a distance of 3.5 Mpc. Its morphology is peculiar; the 
outer isophotes resemble a dwarf elliptical, but over time NGC 5253 has been classifi ed as a spiral, an elliptical, an SO, 
an irregular, and most recently, as an amorphous galaxy (Caldwell & Phillips 1989). A blue starburst dominates the 
central region, with a dust lane bisect i ng the main body along the minor axis. The central starburst comprises at least 
six SSCs, identified bv ICalzetti et al.l (|1997ft . who publis hed HST multiband HST/WF PC2 images of the galaxy with 
a resolution of O '/ l. Th e H I distribution is published bv iKobulnickv fc Skillmanl ([20081 ). with a beam size of 9"0x7"6. 

ICalzetti et al.l (|1997l ) have measured the cluster ages by comparing colors and equivalent width of hydrogen recom- 
bination lines with SB99 predictions. The reddest cluster, N GC 5253-5, has an age of < 2.5 Myr, and dominates the 
infrared spectral energy distribution ([Vanzi fc Sau vagc 2004) ; its visual exti nction Ay is uncertain but could be as 
large as 35 mag (|Calzetti et al.lH997l ). although is probably around 7-8 mag (| Vanzi fc Sauvagell2004f ). The brightest 
and bluest cluster, NGC 5253-4, is also qui te young, ^2.5 Myr, but t he remaining SSCs are o lder, ~10-50Myr. Stellar 
masses of the SSCs have been inferred bv ICalzetti et al.l ([1997) and IVanzi fc Sauvagd (|2004l ). through comparison of 
the observed broadband luminosities to SB99 predictions, given the age of the cluster. Masses range from ~ 10 4 M 
(NGC 5253-4) to 1O 6 M (NGC 5253-5), and radii from 1.6 to 3.5 pc, as measured from HST images. 
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Fig. 10. — Same as Figure [4] but with conservative estimates for the total gas column density inferred from SFRs. (See the electronic 
edition of the Journal for a color version of this figure). 

NGC1140 

NGC1140 is an am orphous, irregular galaxy, and, like NGC5253, has been reclassified over the course of time 
(jHunter et al.l Il994bh . Optically, it is dominated by a supergiant H II regio n encompassing ~ j -0 4 OB stars, far 
exceeding the stellar conten t of the giant H II region, 30 Doradus, in the LMC (|Hunter et al.lll994bl) . The H II region 
is powered by several SSCs (jHunter et al.lll994bt Ide Griis et al.ll2004f). situated in a vertical strip, about 10" in length. 
Hi gh- resolution AST/ W FPC2 images exist for this galaxy (jHunter et all 1199 4b). and H I maps have been published 
bv lHunter et a l. (1994a). The resolution of the H I map (16"x22") is insufficient to distinguish the clusters; given this 
resolution, the peak Shi is certainly underestimated. 

Stellar masses have been determined for the SSCs in NGC1140 by Ide Griis et~aH ((2001 . using a minimization 
technique which simultaneously estimates stellar ages and masses, metallicities and extinction using broadband fluxes. 
We prefer these masses to those measured by a virial technique (jMoll et all I2009J) . as the latter c an be an order of 
magnitude larger, perhaps due to non-virial line widths. Sizes were measured bv [Moll et al.l (j2009l ) for the brightest 
clusters 1 and 6, and we have used these, after correcting for the different distance scale. 

STELLAR MASSES FOR THE LOW-RESOLUTION SAMPLE 

To compute stellar masses for the low-resolution sample of 16 objects, we acquired IRAC 3.6 and 4.5 /jm data from 
the Spitzer archive for all galaxies in the sample (3.6 /jm images are unavailable for two objects). Starting from the 
Basic Calibrated Data, we co add frames using MOPEX , the image mosaicing and source-extraction package provided 
by the Spitzer Science center (jMakovoz fe; Marleaull2005T ) . Pixels flagged by masks are ignored. Additional inconsistent 
pixel values are removed by means of the MOPEX outlier rejection algorithms, in particular the dual-outlier technique, 
together with the multiframe algorithm. We correct the frames for geometrical distortion and then project them onto a 
fiducial coordinate system with pixel sizes of l'/2, roughly equivalent to the original pixels. Standard linear interpolation 
is used for the mosaics. The noise levels in our MOPEX IRAC mosaics are comparable to or lower than those in the 
pipeline products. We then perform aperture photometry on the IRAC images with the IRAF photometry package 
apphot and applying appropriate unit conversion to compute integrated fluxes. The background level is determined 
by averaging several adjacent empty sky regions. Fluxes are computed with a curve-of-growth analysis at radii where 
the growth curve becomes asymptotically flat. From this analysis, we determined the half-light (effective) diameter 
and the size of the object as the point at which the growth curve achieves flatness. Our values for IRAC sizes are 
on average 1.24 times the geometric means of the optical dimensions (as reported in NED). Although the standard 
deviation is large (0.48), this could indicate that there is an evolved stellar population in the extended regions of the 
galaxies that is not seen at optical wavelengths. Alternatively, it could merely be an effect of surface brightness since 
the optical diameters are isophotal and, for a given optical surface brightness, the IR images could be deeper. In 
addition, the two methods of measuring sizes also differ (isophotal in the optical, and photometric in the IR). In fact, 
the approximation of a circular virtual aperture could contribute t o the larg e r size s measured with IRAC. 

From the IRAC photometry, we derive stellar masses following iLee et all (j2006l ) by inferring A-band luminosities 
from IRAC [4.5] total magnitudes (with a color correction), and a color-dependent (B — A") mass-to-light ratio. To 
increase the reliability of this procedure, we include and calibrate with a similar procedure the IRAC [3.6] magnitudes, 
and, where available, also incorporate A-band magnitudes from 2MASS. The average of these three values (one K 
band, and two indirect A-band estimates) are used to calculate the B — A color for the M/L ratio, and the A-band 
luminosity. These magnitudes and colors are reported in Table |3J 

As previously mentioned, inferring stellar masses photometrically can be problematic for some of the BCDs in our 
sample. Hot dust, together with free-free nebular continuum or a high equivalent-width Bra line, can contaminate 
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Fig. 11. — Same as Figure [5] but with upper limits on the total gas surface density inferred from either CO fluxes or SFRs. (See the 
electronic edition of the Journal for a color version of this figure). 
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Fig. 12. — Same as Figure [6] but with lines to connect lower and upper limits on Shi, mod f° r the BR model. Upper limits on Sjjj mo< j 
are derived assuming an upper limit on S gaa . (See the electronic edition of the Journal for a color version of this figure). 



the broadband fluxes from 2 to 5 /xm (|Hunt et al.H2fj0llSinith fc Hancockl [20091 ). For this reason K, [3.6], and [4.5] 
magnitudes can potentially be poor indicators of stellar mass. An extreme case is SBS 0335— 052 E, one of the lowest 
metallicity objects in the sample, where 5 0% of the K band emission is gas, and 13% is dust. Only 37% of the 
2 fim emission is stellar ([Hunt et al.ll200l] ). At 3.8 /im (ground-based L band), the situation is even worse, with 
stars comprising only 6% of the emission. Hence, to mitigate the potential overestimate of the stellar mass from 
contaminated red colors, the minimum (bluest) colors were used to infer the mass-to-light ratio (because of its B — K 
dependence), and the fC-band luminosi ty. A compariso n between the mass-metallicity relation obtained with our 
inferred stellar masses and the sample in lLee et al.1 (|2006f ) suggests that the use of the minimum stellar masses (bluest 
colors) is strongly advocated (with an error of a factor 2-3). 

CONSTRAINTS ON THE TOTAL GAS COLUMN DENSITY 

The analysis presented in the main text is entirely based on lower limits on the total gas surface density, because of 
the impossibility to reliably establish the H2 abundances from the available CO observations. However, we can impose 
conservative upper limits on S gas using indirect ways to quantify Sh 2 ■ 

For the high-resolution sample, we derive S, , s from the molec ular gas column density as inferred by means of the 
SFRs, assuming a depletion time idopi ~ 2 Gyr (jBigiel et al.ll2008l ) for molecular gas. Formally, this should correspond 
to an upper-limit on Sjj 2 , mainly because we use SFRs integrated on scales which are much greater than the individual 
associations we are studying. However, since the H I surface density is not precisely known, we cannot regard S gas 
as real upper limits, although we argue that they likely are. Using these limits, we can explicitly show that the 
disagreement between the extrapolation of the BR model and the observations presented in the Section |4. 1 . 31 cannot 
simply be explained with high gas column densities. This is shown in Figure II 01 where we present once again both 
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the models and data, adding conservative estimates for S gas , connected with a dotted line. The fact that the derived 
S gas ~ 2 - 3 x 10 3 are not enough to account for the observed discrepancy confirm the results inferred using lower 
limits only (Figure |4]). 

Similarly, for the low resolution sample, we set conservative upper limits on E gas assuming a CO-to-H 2 conversion 
factor X = 11 x 10 21 cm -2 K" 1 km" 1 s dLerov et al.ll2009h . Being derived from one the highest X published to date, 
the inferred £h 2 are likely to be truly upper limits on the intrinsic H2. However, whenever these are smaller than 
Sh 2 as obtained from the SFRs combined with a depletion time idcpi ~ 2 Gyr, we assume conservatively the latter 
values. This may be warranted since some galaxies in our sample are at even lower metallicity than the one assumed 
in CO-to-H 2 conversion factor used here. These upper limits on £h 2 are shown in Figure [TT1 and can be used in turn 
to set upper limits on Shi, mod for the BR model (see Figure |6] and Section 14.2. 2|) , as in Figure [12] Because these 
upper limits on Shi, mod exceed significantly the model expectations at low metallicity, we infer that some caution is 
advisable when extra polating local empirical sta r formation laws to high redshift, in dwarf galaxies or in the outskirt 
of spiral galaxies f see iFumagalli fc~Ga vazzi 20081). 
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